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ABSTRACT 


A simple analysis has been developed to describe the transition of 
hybrid rocket combustion from the diffusion-limited region at high pres- 
sure to the kinetic-limited domain at low pressure. Good qualitative 
agreement has been obtained with measured regression rate behavior in 
those operating regimes where surface effects such as excessive melting 
or charring do not occur. 
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I . INTRODUCTION 


Interest in the development of the hybrid rocket as an important 
propulsion device has naturally led to a consideration of its behavior 
over a wide range of operating chamber pressures. Early studies of 
hybrid combustion concentrated on the high-pressure regime in which the 
regression rate of the vaporizing surface is controlled by turbulent 
diffusion in the boundary layer, to the exclusion of any chemical kinetic 
effects . 1-6 The theory shows that in diffusion-limited combustion the 
regression rate is independent of pressure and depends only upon the mass 
flux through the port. 

As the operating pressure is decreased, however, a pressure thresh- 
old is reached below which the chemical reaction time becomes significant 
compared to the mass diffusion time. In this regime of operation the 
regression rate becomes pressure-dependent at a fixed mass flux. Thus, 
the next step in the orderly development of hybrid combustion theory was 
the description of the behavior in the pressure-sensitive regime, which 
was carried out under a previous contract . 7 It was found that a rela- 
tively simple analytical model based on classical turbulent flame theory 
exhibited good agreement with the observed regression rate/pressure 
dependence at a single fixed oxidizer mass flux. 

The current program was designed to supplement the previous investi- 
gation by obtaining data at two other oxidizer mass flux values for 
comparison with analytical predictions and to investigate the growth and 
decay of axial traveling waves produced by injecting pressure pulses 
from an explosive device. Only the regression measurements were carried 
out before cancellation of the contract after a decision was made by NASA 
to terminate their hybrid combustion research contracts. The data obtained 
corroborated the agreement between theory and experiment noted earlier. 
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II. THEORETICAL STUDIES 


Theoretical studies carried out up to the time of the contract 
cancellation were concerned with improving the predictions of the model 
under conditions of varying mass flux. The initial analysis 7 gave good 
agreement with the regression rate/pressure behavior at constant oxidizer 
mass flux, but poor agreement with the regression rate/mass flux behavior 
at constant pressure. The problem has now been corrected by a reconsid- 
eration of the ratio of chemical reaction time to turbulent diffusion 
time in the boundary layer. 

The regression rate r is given by the equation 2 



where r is the regression rate in the diffusion limit, t x is the charac- 
teristic chemical reaction time, and T is the characteristic turbulent 
diffusion time. Note that as t 1 /T-0, r/r -*1; i.e., when the reaction time 
is small compared to the diffusion fime, the diffusion-controlled regres- 
sion rate is recovered. The other limit, corresponding to t 1 /T— 00 , is 
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In this limit, termed the kinetically controlled regime, the flame behaves 
as a premixed flame for which 
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In the original formulation 7 the quantity i x l i 2 appeared as a multi- 
plier on the right-hand side of the equation. The distance was the 
distance that the flame front propagates into unburned gas at the 
(kinetic) flame speed while i 2 was the characteristic scale of turbulence. 
In this formulation t^T was found to be proportional to the square root 
of the expression given in Eq. 4; this gave a dependence on mass flux G 
that did not agree with experimentally observed behavior. Much better 
qualitative agreement with the trends shown by the data is now obtained, 
as will be discussed in the next section. 


To employ Eq. 4 it is necessary to choose a value for the empirical 
constant C. The empiricism enters the problem through the kinetics 
formulation in which the functional dependence of the reaction rate, but 
not its absolute value, is chosen. For the present calculations the 
constant was evaluated by assuming that r/r = 0.75 when p = 150 psia and 
= 0. 05/in. 2 -sec. These numerical choices are supported by the available 
experimental data that are given in the next section. 

Calculations of weight loss as a function of chamber pressure for a 
5-sec time interval are shown in Figs. 1 through 3. Weight loss is shown 
rather than regression rate because it is the measured quantity in the 
experiments. A print out of the computer code used for these calculations 
is given in the Appendix. Because the physical properties (density, heat 
of vaporization, flame temperature, etc.) of PU (polyurethane) and PBAN 
(polybutadiene-acrylic-nitrile) are nearly the same, 7 the computed weight 
losses for the two binder systems differed by no more than 2 °jo over the 
ranges of pressure and oxidizer mass flux considered. Therefore, there 
is no differentiation between the two binder systems delineated in the 
figures. Oxygen was considered as the oxidizer in all of the calculations. 
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WEIGHT LOSS AT 5 SEC — grams 
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FIGURE 1 THEORETICAL WEIGHT LOSS AS A FUNCTION OF CHAMBER PRESSURE 
(G° = 0.05 lb/in. 2 -sec) 
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WEIGHT LOSS AT 5 SEC — grams 
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THEORETICAL WEIGHT LOSS AS A FUNCTION OF CHAMBER PRESSURE 
(G° = 0.10 Ib/in. 2 — sec) 
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FIGURE 3 THEORETICAL WEIGHT LOSS AS A FUNCTION OF CHAMBER PRESSURE 
<G° = 0.15 lb/in. 2 -sec) 
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Figures 1 through 3 show that the dependence of weight loss on 
chamber pressure at low pressures becomes more marked as the oxidizer 
mass flux increases. Increasing the metal loading at a fixed oxidizer 
mass flux also amplifies the effect, in this case because of the depen- 
dence of radiative heat transfer on pressure. 

Figures 4 through 6 are cross-plots of the results shown in Figs. 

1 through 3 to show the dependence of weight loss on oxidizer mass flux 
as a function of chamber pressure. At low chamber pressures, in the 
kinetically controlled regime, the weight loss is nearly independent of 
oxidizer mass flux. As the pressure increases the slope of the curve 
increases, with the proportionality approaching something less than the 
0.8 power. Here the diffusion-controlled regime, in which regression 
rate is proportional to the 0.8 power of the total mass flux, is being 
approached. 

The trends predicted by the theory are, of course, those that were 
built in through the original formulation of the physical problem. The 
problem was formulated on the basis of past results obtained at SRI 7 and 
elsewhere. 8-10 A comparison with more recent data obtained under this 
contract is given in the next section. 
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FIGURE 4 THEORETICAL WEIGHT LOSS OF A PURE BINDER HYBRID GRAIN AS 
A FUNCTION OF OXIDIZER MASS FLUX 
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FIGURE 5 THEORETICAL WEIGHT LOSS OF A HYBRID GRAIN CONTAINING 20 PERCENT 
ALUMINUM AS A FUNCTION OF OXIDIZER MASS FLUX 
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FIGURE 6 THEORETICAL WEIGHT LOSS OF A HYBRID GRAIN CONTAINING 40 PERCENT 
ALUMINUM AS A FUNCTION OF OXIDIZER MASS FLUX 




III. EXPERIMENTAL STUDIES 


The experimental studies were designed to provide data that would 
assist in the development of a theoretical model that could be used to 
predict the regression rate behavior and combustion instability limits 
for hybrid rocket motors operating in the low-pressure regime. The 
experiments originally proposed for this program included four types of 
tests: 

1. Hybrid motor tests for determining the mean regression rate of 
various propellant systems in the pressure-sensitive regime over 
a range of oxidizer mass flux. 

2. Hybrid motor tests in which axial instability would be induced 
by pressure pulses. The growth and decay of the pressure wave 
would be studied under varying operating conditions. 

3. Slab burner tests to obtain Schlieren photographs of the rela- 
tive thickness and position of the flame zone under the various 
operating conditions chosen. 

4. Atmospheric combustion tests in the combustion simulator to 
determine the effect of heat release distribution on the heat 
transfer to the wall. 


Because of the premature termination of the program only the first phase 
of the experimental studies was compileted. 


A. Apparatus for Regression Rate Measurements 

The hybrid motor tests for determination of the mean regression 
rate utilized the existing flow facility that is shown schematically 
in Fig. 7. This facility was originally designed for liquid rocket 
engine tests with cryogenic oxidizers and was therefore constructed of 
304 stainless steel with Teflon* and/or Kel-F^ seats and seals. Sonic 
chokes operated above their critical pressure ratio were used to maintain 
a constant oxidizer mass flow throughout each test. These chokes were 
calibrated against a flow meter in the preceding hybrid combustion insta- 
bility program (Contract No. NAS 1-7310). 7 

*Trademark, E. I. du Pont de Nemours and Co. 
tTrademark, 3M Company 
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The motor configuration used was the same as in the prior program. 
The 2-l/2-inch-diameter motor consisted of a 12-inch flow-straightening 
section in which the ignition system was housed, a 12-inch tubular grain 
with a 1-inch-diameter internal perforation, and a water-cooled nozzle 
assembly. 

The grain is ignited by preflowing a small amount of oxygen and 
propane into the chamber and igniting it with a spark plug. The ignition 
system is preset so that it operates near the lower flammability limit 
for the oxygen-propane mixture, thus yielding an oxygen-rich high- 
temperature source. The time for the flame to spread from the head end 
of the grain to the nozzle end was determined to be approximately 1 sec. 
Therefore, the main flow valve was delayed for that length of time. The 
weight of fuel typically consumed during this ignition process was de- 
termined to be 0.5 g which is about 1 °jo of the amount of fuel that would 
be consumed in a 5-sec test at a regression rate of 0.01 in. /sec. 

B. Regression Rate Measurements 

Weight loss data were obtained for test durations of both 5 and 10 
seconds. The grains were dimensionally measured before and after each 
test. The data were obtained for PU and PBAN binders containing 0, 20, 
and 40 percent aluminum for chamber pressures ranging from 30 to 250 psia 
and for oxidizer mass fluxes of 0. OS', 0.10, and 0.15 lb/in. 2 -sec. The 
oxidizers used for all of the tests was gaseous oxygen. The data from 
these tests are presented in Table I. PU and PBAN were chosen for com- 
parative binder tests because their surface behavior during regression 
is quite different. PU sublimes cleanly at the regressing surface after 
a melt is formed, while the regression of PBAN consists of the stripping 
of hydrogen from the surface layer, leaving a carbonaceous char that sub- 
sequently breaks away in particulate form. 

Since these two systems behave differently the weight loss data will 
be subsequently discussed separately for each system; however, a brief 
discussion of the data handling procedure is in order. 
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Table I 


EXPERIMENTAL DATA SUMMARY 


Grain 

Number 

Binder 

Type 

Aluminum 

($) 

Chamber 

Pressure 0 

(psia) (lb/in. 2 -sec) 

Weight 

Loss 

(grams) 

Run 

Duration 

(sec) 

508 

PU 

0 

34 

0.05 

45 

5. 1 

509 

PU 

0 

33 

0.05 

86 

10.0 

510 

PU 

0 

77 

0.05 

58 

5.0 

511 

PU 

0 

78 

0.05 

107 

9.6 

512 

PU 

0 

165 

0. 05 

77 

5.0 

513 

PU 

0 

170 

0.05 

147 

9.4 

514 

PU 

0 

260 

0.05 

79 

4.0 

515 

PU 

0 

253 

0.05 

156 

9.2 

527 

PU 

0 

75 

0.10 

59 

5.0 

528 

PU 

0 

260 

0. 10 

97 

5. 0 

517 

PU 

0 

35 

0. 15 

60 

5.0 

518 

PU 

0 

35 

0.15 

108 

9.6 

519 

PU 

0 

80 

0. 15 

67 

5.0 

521 

PU 

0 

74 

0.15 

124 

9.6 

523 

PU 

0 

160 

0. 15 

87 

4.9 

524 

PU 

0 

130 

0.15 

163 

9.6 

525 

PU 

0 

225 

0.15 

92 

4.9 

526 

PU 

0 

230 

0.15 

202 

9.6 

556 

PU 

20 

34 

0.05 

49 

4.7 

557 

PU 

20 

35 

0.05 


9.5 

558 

PU 

20 

72 

0.05 

60 

5.0 

559 

PU 

20 

81 

0.05 

133 

9.8 

560 

PU 

20 

130 ^ 

0.05 

82 

5.0 

561 

PU 

20 

138 

0.05 

174 

9.6 

543 

PU 

20 

36 

0. 10 

59 

5.0 

544 

PU 

20 

35 

0.10 

128 

9.8 

545 

PU 

20 

71 

0. 10 

84 

5.0 

546 

PU 

20 

87 

0. 10 

186 

9.8 

547 

PU 

20 

156 

0. 10 

130 

5. 0 

548 

PU 

20 

163 


252 

9.8 

549 

PU 

20 

32 

0. 15 

65 

5.2 

550 

PU 

20 

35 

0. 15 

128 

9.7 

551 

PU 

20 

62 

0. 15 

82 

5.0 

552 

PU 

20 

83 

0. 15 

196 

9.8 

553 

PU 

20 

157 

0. 15 

139 

5. 1 

554 

PU 

20 

169 

0. 15 

277 

9.5 


14 




Table I (contd) 
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Table I (concld) 
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Typical weight loss data as a function of time are presented in 
Fig. 8. The mean chamber pressure at which the fuel is burnt is denoted 
beside each data point. Since minor variations in test duration did 
occur, the data were standardized at 5 sec. This relatively short time 
was chosen to minimize the effect of the change in oxidizer mass flux 
as the internal perforation diameter increased. The data were then 
plotted as a function of chamber pressure and oxidizer mass flux for the 
various aluminum loadings made. From the weight loss versus chamber 
pressure curves, constant pressure lines were constructed for the weight 
loss verus oxidizer mass flux curves. The chamber pressure values used 
are mean values. For either pure PBAN or PU the observed pressure devi- 
ation about the mean was very minor. However, in the systems that con- 
tained aluminum the deviations were much more severe because of the coating 
of the nozzle throat with aluminum oxide. The most severe cases were those 
in which the chamber pressure was high and/or the oxidizer mass flux was 
low, both cases necessitate small nozzle throat sizes. With this varia- 
tion a precise mean value of chamber pressure for repetitive tests were 
difficult to obtain. 

Under the current program the first data to be taken were the weight 
loss data at the low value of oxidizer mass flux. These data are pre- 
sented in Fig. 9. Upon examining the data, it can be seen that the PBAN 
results fall off for chamber pressures above 150. One possible explanation 
is that this decline occurs when the velocity is reduced below the level 
required to efficiently remove the char layer that is present. This idea 
is supported by the results shown in Fig. 10 where the oxidizer mass flux 
is 0.10 lb/ in. 2 ~sec. Here the velocity is large enough to stabilize the 
regression rate for pressures over 150. In Fig. 11 the oxidizer mass flux 
is further increased and the regression rate becomes a function of chamber 
pressure throughout the range of test conditions. 

The dependency of weight loss of the PBAN system on oxidizer mass flux 
is depicted in Figs. 12 through 14. At lower values of chamber pressure 
the regression rate is relatively independent of mass flux. However, as 
the chamber pressure is increased the transition from a kinetically 
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WEIGHT LOSS — grams 



FIGURE 8 TYPICAL WEIGHT LOSS OF A HYBRID GRAIN AS A FUNCTION OF TIME 
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WEIGHT LOSS AT 5 SEC — grams 



FIGURE 9 EXPERIMENTAL WEIGHT LOSS OF HYBRID GRAINS AS A FUNCTION OF 
CHAMBER PRESSURE (G° = 0.05 lb/in. 2 -sec) 
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WEIGHT LOSS AT 5 SEC — grams 



FIGURE 10 EXPERIMENTAL WEIGHT LOSS OF HYBRID GRAINS AS A FUNCTION OF 
CHAMBER PRESSURE (G° = 0.10 lb/in. 2 -sec) 
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WEIGHT LOSS AT 5 SEC — grams 



FIGURE 11 EXPERIMENTAL WEIGHT LOSS OF HYBRID GRAINS AS A FUNCTION OF 
CHAMBER PRESSURE (G° = 0.15 lb/in. 2 -sec) 
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WEIGHT LOSS AT 5 SEC — grams 



FIGURE 13 EXPERIMENTAL WEIGHT LOSS OF A HYBRID GRAIN CONTAINING 20 
PERCENT ALUMINUM AS A FUNCTION OF OXIDIZER MASS FLUX 




200 



FIGURE 14 EXPERIMENTAL WEIGHT LOSS OF A HYBRID GRAIN CONTAINING 40 
PERCENT ALUMINUM AS A FUNCTION OF OXIDIZER MASS FLUX 
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controlled regression rate to a diffusion-limited process is apparent. 

As mentioned previously the PU system exhibits a different combustion 
mechanism than the PBAN system. As can be seen from the weight loss versus 
chamber pressure data which are also shown in Figs. 9 through 11, the PU 
system's regression rate is dependent upon chamber pressure for all of the 
values of oxidizer mass flux investigated. Since this system melts rather 
than chars, the velocity effect at constant oxidizer mass flux is absent. 
However, the influence of velocity is noticeable at low chamber pressures 
and high oxidizer mass flux. This is thought to be an erosion effect which 
is diminished by the addition of the aluminum. The aluminum not only in- 
creases the viscosity of the melted surface layer but also reduces the 
thermal profile thickness by increasing the regression rate. 

C. Comparison with Calculated Weight Loss Behavior 

A comparison between Figs. 1 through 3 and Figs. 9 through 11 shows 
that the qualitative behavior of weight loss as a function of chamber 
pressure is in best agreement with the calculated behavior at the highest 
oxidizer mass flux. At the lowest value of mass flux the PBAN curves 
tend to droop at high pressure, presumably because the velocity is too 
low to remove the char layer, and the PU curves tend to rise at low 
pressure, presumably because of erosive effects on the melted binder surface 
at low pressure. This behavior is also reflected in a comparison of Figs. 

4 through 6 with Figs. 12 through 14. 

In all cases the measured weight loss is larger than the predicted 
value. The reason for this discrepancy was not known at the time work was 
halted on the contract. 
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IV. CONCLUDING REMARKS 


A simple analysis to describe the transition of hybrid rocket com- 
bustion from the diffusion-limited region at high pressure to the 
kinetic limited domain at low pressure has been developed. Good quali- 
tative agreement has been obtained with measured regression rate behavior 
in those operating regimes where surface effects such as excessive melting 
or charring do not occur. 

It remains to be seen whether the simple model can describe insta- 
bility behavior in the pressure-sensitive regime. Future investigations 
of hybrid combustion phenomena should be directed toward an understanding 
of this regime of operation. 
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Appendix 

COMPUTER CODE FOR CALCULATION OF HYBRID REGRESSION RATES 
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PROGRAM HYBRID(INPUT, OUTPUT, TAPE5*INPUT*TAPE6=0UTPUT> 001 

000003 PEAL l»K»*1 *k2*k3.K4,k5*W8,K9*M,mO.M] ,m2,h3.N( 100) .NN.mH 00? 

000003 DIMENSION 0( 100) *OATE < 12) «FUEL (12) *0XIn(12> no"* 

C THIS PROGRAM IS BASED UPON MEAT-TRANSFER limited HYBRID theory, 004 

C IT CALCULATES THE REGRESSION BEHAVIOR OF A CIRCULAR TUBE FUEL 005 

C GRAIN. THE INPUTS WHICH ARE NECESSARY FOR OPERATION OF THE 006 

C PROGRAM 4fi£ 0EF1NE0 AS FOLLOHSI * 007 

C 00 IS INITIAL INNER DIAMETER (IN.) 000 

C L IS LENGTH (IN.) OOP 

C U2 IS OUTER DIAMETER (IN.) 0)0 

C K IS HASS FRACTION OF NON-VOLATILE SURFACE MATERIAL "11 

C K 1 IS MASS FRACTION OF VOLATILE SURFACE MATERIAL WHICH FORMS 012 

C PARTICULATE COMBUSTION PRODUCTS 013 

C 0 ? IS TOTAL OENSITY OF FUEL grain |LB/IN<**3) 014 

C *3 IS EFFECTIVE HEAT of GASIFICATION (BTU/LB1 015 

C MO IS OXIOIZER FLO* RATE (LB/SEC) 016 

C (IF MO VARIES* IT HAY EITHER BE REA0 AS DATA OR 017 

C DESCRIBED BY AN EQUATION.) Olo 

C PO IS CHAMBER PRESSURE (PSU) (SEE THE ABOVE) 119 

C j IS CASE NUMBER 0?0 

C J«1 IF PARTICLES IN THE GRAIN PRODUCE PARTICULATE 021 

C COMBUSTION PRODUCTS 0?2 

C J«2 IF PARTICLES IN THE GRAIN PRODUCE GASEOUS 023 

C COMBUSTION PRODUCTS 024 

C J«3 IF A COMPLETELY VOLATILE GRAIN PROOUCES PARTICULATE 025 

C COMBUSTION PRODUCTS 0?6 

C UO IS OXIOIZER TEMP. AT INLET (OEG RaNKINE) 0?7 

c us is combustion temp, at an o/f ratio which is 3/4 of the o?a 

c stoichiometric value AT 5 ATM. PRESSURE (DEG RANAINE) 0?q 

c C7 is gas-phase radiation constant 030 

C K* IS MASS OF OXIOIZER CONSUMED per unit MASS OF NON-VOLaTiLE 031 

c surface material in producing particulate combustion products 032 

C K5 is MASS OF OXIOIZER CONSUMEO PER UNIT MASS OF VOLATILE 033 

c surface material in producing particulate combustion products 034 

C C8 IS MASS OF PARTICULATE COMBUSTION PRODUCTS FORME0 pen unit 035 

C MASS OF NON-VOLaTILE SURFACE MATERIAL 036 

C E2 IS EMISSIVITY ( ABS0RPT I V ITT I OF FUEL SURFACE 037 

C B IS HASS TRANSFER NUMBER 

C K8 IS OENSITY OF PARTICULATE COMBUSTION PRODUCTS (LB/IN*»3) 039 

C C2 IS BOUNDARY LAYER DENSITY CORRECTION FACTOR O 4 O 

C S IS OXIOIZER VISCOSITY AT TEMP. Uo (LB/IN-SEC) OA) 

C T 1 IS TIME INCREMENT FOR CALCULATIONS (SEC> 042 

C T 2 IS QUIT TIME (SEC) 043 

C T 3 IS TIME INCREMENT FOR PRINTOUT (SEC) 044 

C XI IS OISTANCE INCREMENT FOR CALCULATIONS «X/L) 0a5 

C C9 IS ThE RAOIATION TEMPERATURE CONSTANT 

C R6 IS THE PARTICLE RAOIUS IN MICROnS 

000003 RE AO (5*1000) (OA TE ( I > * I»1 • 12 > 046 

000015 1000 FORMAT ( 12A6) 047 

000015 10 READ (5* 1010) DO*L*02*K*K1*K2*K3*K4«K5*C8*C9*T1*T2*T3* 042 

1 M0*P0*U0*U5*E2«B*K8.R6«QF 049 

000077 IF < 00*EQ* 0 • 0 ) GO TO 850 050 

000100 1010 FORMAT (7FlO*4) OSl 

000100 REAO (5*10201 J 05 ? 

000106 1020 FORMAT (II) 053 

000106 REAO (5. 1000) (FUEL < I ) * I»1 . 12 ) 054 

000120 REAO (5 * JOOO I (0X10(1) *1*1*12) 055 
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000132 

E3*0.0 


056 

000133 

U2*L5*(1.0*3.33E-*»<PO-75.0> > 


057 

OOOUO 

C7*0,003 


059 

0001*1 

CP=0.* 



0001*3 

NN=2.0 



0001** 

£F*27000.0 



0001*6 

RR*0.75 



0001*7 

TEH*2.0 



000151 

3000 T£X*2.0*T£R-RR**2*TER«»2.2.0*EXP t-TER> 



000162 

IF ( TEX.LT.2. Ol ) GO TO 3020 



000166 

TERMER. 0.01 



000167 

60 TO 3000 



000167 

3010 TERMER-o.Ol 



000171 

60 TO 3000 



000172 

3020 IFHEX.LT. 1.09)60 TO 3010 



000175 

■RITE(6,2070)TER 



000202 

HR"0 . 0*33 



000203 

AR*0 .6659 



000205 

PR*1SO.O 



000206 

U"R»1**0.0 



000210 

UFH«U5*(1.0«3.33£-**(PR-7S.0> ) 



000215 

£l"3.307E-15 


0*0 

000216 

C6*(3.0*16.387£12|/(*.0* K 8*R6**3> 



000222 

Jl"l 


063 

000223 

Cl *0.036 


06* 

00022S 

C2*1.0 


065 

000226 

S*3.00E-6 


066 

000230 

C3» ( S/L) **0 .2 


067 

00023* 

C5»3. 1*159 


068 

000235 

C*"C5/*.o 


069 

000237 

Xl*0.02 


070 

0002*0 

I*»INT(0.1/X1»0.5) 



0002** 

**0.0 


07? 

0002** 

U1*C9*U2 


073 

0002*7 

"RITE (6.8888) 


07* 

000252 

8888 FORMAT (1H1) 


075 

000252 

WRITE (6.2000) c DATE 1 1 ) . 1*1. 12> 


076 

00026* 

2000 FOHMATIIH .50X.12A6/) 


077 

00026* 

WHITE (6,2010) (FUEL (I)»I*1.12) « 1 OX ID < I ) « I = 1 « 1 2 ) 


078 

00030* 

2010 format (5oH REGRESSION BEHAVIOR OF A CIRCULAR TUBE FUEL GRAIN// 

079 


16H FUEL*. 12A6/10H 0XI0I ZER* » 12A6/ ) 


080 

00030* 

WRITE (6,2020) 00 ,L« K2 ,k3 ,B,PO ,mO ,uO ,u 2 ,Ul ,C7 ,f2 , K . K 1 . K* »K5 . C8 

OrI 

000352 

WRITE (6«2o21 ) t1,x1,t3,T2,J 


08? 

000370 

2020 FORMAT U9H INITIAL CONOITIONs/5>15H oiameter 

« ,E 1* ,5/ 

083 


1 15H LENGTH *«El*.5'l5H DENSITY 

= 'EH.5^ 

08* 


2 15H EFF. H-5UB-V *.El*.5/lSH BEE 

*.E1*,5? 

085 


3 ISm PRESSURE *,E1*.5/15h OXIq. FLOW 

**E 1* .5/ 

086 


* 15H INLET TEHP. *.E 1* .5/ 15H FLAmE TEMP, 

s,E 1* ,5/ 

087 


5 15H RAO. TEMP. *.El*.5/15M GAS CONST. 

».El*.5? 

088 


6 15M WALL ABSORPT .*,El*.5/lSH MgT*L PCT. 

*,E 1* .5/ 

089 


7 15H PART. PCT. «,El*.5/l5H ZET* 

»,El* .5/ 

090 


8 15H ZETA-1 «,E1*.S/1SH LAM 8 DA 

*,E 1*,5) 

091 

000370 

2021 FORMAT <15R OEL-T »,El*.5/l5M OEl-*/L 

**El*.5? 



1 1SH PRINT TIME *«El*.5/15M 6UIT TIME 

2 15M CASE *.!*) 

",E 1* ,5/ 


000370 

WRITE (6,8888) 


092 

00037* 

WR1TE(6,2000) (DATE! I) ,1*1 ,12) 


093 

000*06 

U6*C 1*C2*C3 # <8 ) **0. 23 


09* 


29 



005 


000*1* 


T**0.0 

000*15 


T5*0.0 

000*16 


T«0.Q 

000*20 

155 

*»*l 

000*21 


1*1 

000*22 


11*1* 

000*2* 


M*M0 

000*26 

165 

IFU.GT.O.OIGO TO 175 

000*31 


0(1) s 00 

000*33 

175 

R3*(X*L)/(S. 0*0(1) ) 

000*37 


IF ( 83 .UE • 1 • 0 > GO TO 190 

000**2 


IF <83. GT. 5.0)00 TO 200 

000**5 


IFI83.GT.1.0IGO TO 210 

000**7 

190 

E»0<I)»(1. 0-0.21 •( 83 >»«0.8> 

000*56 


U»U0 

000*60 


GO TO 220 

000*61 

200 

E*0(l)*(1.0-0.21) 

000*6* 


U«U2 

000*66 


GO TO 220 

000*67 

210 

E-D < l J • ( 1 .0-0.21 ) 

000*72 


U*UO«( (u2-U0)/*.0)*(H3-1.O) 

000501 

220 

fl*C**(E>**2 

000503 


Q1*C6*(X)«« (-0.2) •<H/*I ••0.8 

000515 


Q3*C6*01 , 00*X1 , L 

000522 


81*C1/(K2*(1.0-K) 1 

000526 


GO To (255.3*5.3001 .J 

000535 

255 

IFd.GT.llGO TO 280 

0005*1 


M*M0*<1.0-8*8*/( 1.0-8) 1*03 

000550 


Ml=M0*Q3«1.0/<1.0-8) 

000555 


Q5*03 

000556 


N ( 1 1 *C6*C8*8/ (1*0-8)*PO/U*1.0/M*03*(R3)**(-0,8) 

00057* 


GO TO 366 

000575 

280 

Mmm* ( 1 .0-8*8*/ ( 1 .0-8) 1 *G* 

00060* 


m1»m1.Q* *1,0/11.0-8) 

000612 


IF (83 «GT , 1 . 0 1 GO TO 295 

000615 


N(1>«C6*C8*8/<1.0-K)*PO/U*1,0/M*Q5*(R3>**(-0.8> 

000633 


GO TO 366 

00063* 

295 

N(I)-C6«C8*8/(l.o-8)*PO/U*l.o/R»05 

0006*6 


GO TO 366 

000646 

300 

IFd.GT.llGO TO 325 

000652 


M-80*<1.0— Kl— 51«K5)*Q3 

000657 


Ml *MU *03 

000660 


05*03 

000661 


N ( 1 1 *C6*8 1*PQ/U*l , 0/M*Q3* (R3 1 •# (-0 .8 > 

00067* 


GO TO 366 

000675 

325 

M*M* ( 1 .0-81-81*85) *0* 

000702 


Ml*Ml «Q* 

00070* 


IF ( 83.GT . 1 . 0 ) GO TO 3*0 

000707 


Nd)*C6*81*P0'U*l ,0/M*Q5*(R3>** (-0.8) 

000721 


GO TO 366 

000722 

3*0 

N(I)*C6*81*P0/U*1 ,0/M»G5 

000731 


GO TO 366 

000731 

3*5 

IFd.GT.llGO TO 360 

000735 


M*MO«( 1,0/ (1.0-8) 1*03 

0007*1 


M 1 *F 

0007*2 


Nil) *C7*P0 

0007*5 


GO TO 366 

0007*6 

360 

M«M*(1, 0/(1. 0-8) )•«* 


006 

007 

008 

000 
100 
lnl 
10? 
103 
10 * 
105 

1 06 
107 

1 OR 

109 

110 
111 
11 ? 
113 
11 * 
115 
115 
11 ? 
118 
119 
1?0 
>?1 
1?? 
1?? 
1?* 
1?5 


1?8 

1?9 

130 


* -'C 

133 

13* 

135 

136 

137 

138 

139 
1*0 
1*1 
1*? 
1*3 
1 ** 
1*5 
1*6 
1*7 
1*« 
1*9 

150 

151 
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000752 Ml-P 152 

000753 N(I)-C7»P0 153 

000756 366 IF ( jl .EQ.2 ) SO TO 371 Is* 

000760 R2»R1 155 

000761 ua-OaR/ < 1 .0-U8R-AL0G ( 1 0 . 0-R2 ) /EF ) 

000771 TEE»TER» (a/HR»AR/A« (PR/P01 ••NN , IUFR/U2)* # (NN«2.0)» I 1.0-CP»(UFR-UaR 

1 )/0F)/ (1.0-CP*(U2-Ua)/QF)-EXP(-UFR/U2> ) 

001032 R2-R2-I <2.0/TEE)**0.5-(1.0-l .0/TEE-<1 .o-EXP(-TEE> ) >*-0.5) 

001051 SO 70 390 156 

001052 371 IF(J. £0.2)80 70 38* 157 

00105* IF(MI).G7.0.0)G0 70 375 158 

001057 R2»R1 159 

001057 Ua-liaR/ ( 1 . 0-UaR-ALOG ( 10 .0-R2 > /EF ) 

001067 7EE»TER*<P/MR*AR/A«(PR/P0)»«NN»(UFR/u2) , »tNN*2.0)»U.0-CP»(UFR-UaR 

1 )/QF)/(1.0-CP«(u2-oa)/QFI»EXP(-UFR/U2> ) 

001130 R2-R2*< (2.0/TE£)« # 0.5-(1.0-1.0/TEE«< l .o-EXP(-TEE) I >••0.51 

0011*7 GO TO 390 160 

001150 375 IF(I.LT«I*»60 70 382 

001153 IF (I. GT. <6«1*> )G0 TO 382 

001156 IF(I.EQ.Il)GO TO 382 1*3 

001157 GO TO 387 16* 

001157 382 IF ( T.EQ» T5 ) GO TO 590 

001161 GO TO 387 166 

001162 38* IF1I.GT. <2*1*1 1 00 TO 387 

001166 GO TO 375 168 

001166 387 Q2*E 1 # E2/K3 , E3* <U1 ) 170 

001173 H2-Ql*(02/0l*£XP(-a2/Ql) )/( k 2*(1.0-K) ) l7l 

001206 Ua-liaR/ ( 1 . 0-U8R-ALOG ( 1 0 . 0-R2 ) /EF ) 

001215 TEEaTER* (P/MR«AR/A* (PR/Po»»«NN* (UFR/U2) •• (NN*2»0) • 1 1 . O-CP- ( UFR-uaR 

1 )/QF)/(1.0-CP«<U2-Ua>/QF)*EXP<-UFR/U2> ) 

001256 R2«R2*( (2. O/TEE 1«»0.5« (1.0-1. 0/TEE»( 1.0-EXP (-TEE) ) >••0.5) 

001275 390 0*bC5*K2«(1.0-K)»R2*0(I)«X1*L 172 

001305 05-05*0* 173 

001306 IF ( I ,EQ. 1 1 ) GO TO *95 )7* 

001311 D ( 1 ) -0 ( I ) «2.0*R2*T 1 175 

00131* IF (0 ( I ) .67 .02)00 TO 795 1 76 

001320 *05 IF(X.GT. 0.995) GO TO *25 177 

00132* 1-1*1 178 

001325 X-X«X1 179 

001327 GO TO 165 laO 

001327 *25 IFIT.EO. 0.0>GO TO *50 1-1 

001330 m3- (H 1 *M2 ) /2 • Q-mO la2 

00133* a-*«83-Tl 1 fl3 

001336 m2-h1 1r* 

001337 IFIT.EO. T*)GO TO *55 IRS 

0013*2 T-T* T l ** 186 

0013*2 IFIT.gT.T2) GO TO 800 187 

0013*6 GO TO 155 lap 

0013*6 *50 *2*8 1 189 

0013*7 IF (01 . EQ.2 ) GO TO *55 Ion 

001352 Jl-2 l9l 

001353 GO TO 155 192 

001353 *55 aRITE (6.2030) 193 

001357 2030 F0R8ATI/13H AVG. 0EL-R.6X.8HaT. L0SS.6X .SHGAS FLOa.SX. 19* 

l 10HTOTAL FLOa. 7X . 3HG-0 ) 195 

001357 aRITE (6.203S) 

001363 2035 FORPA T ( I OH { IN) . 1 1 X «*H (GM ) »8X .8h (lB/SEC) »6X .8M (LB/SEC ) . 

1 3X.12H(L8/IN2-SEC) ) 
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001363 

001373 

001376 

00U00 

00140* 

0014J1 

001421 

001423 

001425 

001430 


001430 

001432 

001437 

001441 

001441 

001443 


001444 

001446 

00144? 

001455 

001455 

001461 

001461 

001465 


001465 

001471 

001*71 

001474 

001477 

001501 

001522 

001522 

00152? 

001533 

00153* 

001535 

001540 

001542 

001544 

001545 

001556 

001557 

001561 

001562 

001565 

001566 

001570 

00157* 

001600 

00161* 

001617 

001622 

001625 

001626 

001627 

001631 


06»(»/(C44K2»L)« (001«*2)4»0.5 

AVGO* (06-00 1/2.0 

MTL««»453.6 

GSUBO*mO/(C*»(06I«*2) 

UNITE (6«2040> AVGD.uTL.M.Ml.GSUBO 
20*0 FOKFAT (5E14.5//I 
TuT ♦ T 1 
T*«74«T3 

IF ( T .GT . T2 ) GO TO 800 
GO TO 155 

495 IF (T.EQ<T*)GO TO 506 
0(11 s 0 (It <2 .0«R2*T1 
ll»Il<I* 

GO TO 405 

506 IF ( J1 .EQ.2 ) GO TO 510 
GO TO *05 

510 IF ( 1 1 .EQ. I* 1 GO TO 520 
GO TO 535 

520 UNITE (6 <2050 ) T 
2050 FONMAT(/llH TIME/FlO.O/) 

UNITE (6<205S) 

2055 FOHFATI25H TOTAL < 7*. 1 OMCONVECT I VE ) 

UNITE (6.2060) 

2060 FONFAT ( 1 3h 0ISTANCE.5X.8HN EG NATE .6X.8HREG NATE.SX. 

1 5H0EL-H.7X.9MTUBE AREA,5X.9HFL0u ANE A. 5 X < 1 OHEM ISSI VI T Y ) 

UNITE (6<206S ) 

2065 FONFAT ( 1 1 H (X/L) <7X,8HUN/SEC) <6X,8M( IN/SEC) <8X,4H(IN) < 

1 9X,7M(IN«u 2) <7X.7h(IN4*2)/) 

535 DELN«(D(I)-00)/2.0 
TUBEAbC** (0(1) ) **2 
FL0UA-C*»<E) ,# 2 

UNITE (6.20 70 I *<R2 <Nl <OELR< TU8EA.FL0UA.E3 
2070 F0NFATI7E14.5) 

D ( I ) >0 ( I ) <2, 0*N2*T 1 
IF(Dd) .GT.021G0 TO 795 
Il«Il»I4 
GO TO *05 
590 T5*T5* (2.0*T3I 
600 N«U (11/2.0 
N0-N-0.25 
YluS. 1*159/30.0 

Y2»AT AN (1*0/ ( (R»R)/ <RO«R01-1.0>»*0.5) 

E3«0 .0 

J5-1 

633 K9«1.0 

YuYl/100.0 
635 Z>3. 14159/2.0 

J4u 1 

645 Y3u2.0»H«C0S(T) 

IFIY.GT.Y2) GO TO 650 

S3* ( Y3/2. 0- < N0uR0-R«R4 ( SIN ( Y ) >*»2>*»0.5>/SIN(Z> 

650 X3»J4*X1«L 

S4»Y3/SIN(Z) 

IFlY.LT.Y2) GO TO 666 

M9»OiO 

S3«0.0 

666 IF < u5 .EQ. 1 1 GO TO 670 

S5UK9«EXP(-MI-J*I*S3)-K9*E*P<-NII-J4)»<S*-S3) ) «EXP I -N ( I- J4 ) *S4 ) 


196 
1 07 
1 98 
1 99 
200 
201 

212 
2(V3 

264 

205 

206 
207 
2"R 

209 

210 
211 
212 

213 
21* 
2)5 


216 

217 

219 


2 1 9 
220 
221 


22 * 

226 

226 

227 

229 

229 

230 

231 

232 

233 
23* 

235 

236 

237 

238 

239 
2*0 
2*1 
2*2 
2*3 
2 ** 
2*5 
2*6 
2*7 
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001656 GO TO 675 2*8 
001657 670 S5*K9*EXP < -N ( I ♦ J* ) *S3 ) ”K9*EXP (-M I ♦ J* ) * t S*”S3 I > »EXP ( «N ( l ♦ J* ) »S* ) ;>*<> 
00170* 675 S6»<1.0-SS)*<SIN(Z> )«*2*C0SIY> ?50 
00171* Zl*(Y3*Xl*L)/< (X3) **2» < Y3 ) »»2 ) 251 
001721 S7»S6*Z1»Y1*2. 0/3. 1*159 252 
001724 E3*£3«S7/1. 06988 253 
001727 Z«Z-Zl 25* 
001731 J*»J**1 255 
001733 IF < J5.E0. 1 ) GO TO 705 256 
001735 IF(j*.GT.IIGO TO 720 257 
0017*0 IF ( J* «GT . 1 0 ) GO TO 720 258 
0017*2 GO TO 6*5 259 
0017*2 70S IF(d*.GT. (1.0/X1-I))G0 TO 720 260 
001751 IF(J*.GT.10)GO TO 720 2 6 1 
00175* 60 TO 6*5 262 
00175* 720 Y»Y»Y1 261 
001756 IFIY.GT. (3. 1*159/2. 0I1G0 TO 7*5 26* 
001761 GO TO 635 265 
001761 7*5 IF ( JS.EU.21G0 TO 750 266 
001763 J5«2 2*7 
U0176* GO TO 633 2*8 
001765 750 GO TO 387 269 
001766 795 *«ITE (6.2080 • 27n 
001772 2080 FORMATI/37H SPECIFIED EXTERNAL DIAMETER EXCEEDED/) ?T1 
001772 800 GO TO 10 272 
001773 850 STOP 273 
001775 ENU 27* 




33 



REGRESSION BEHAVIOR OF A CIRCULAR TUBE FUEL GRAIN 


FUELwPU 

OXIOI2ER»OXYG£N 
INITIAL CONDITIONS 


DIAMETER 
LENGTH 
DENSITY 
EFF. H-SU8-V 
BEE 

pressure 

OXIO. FLOW 
INLET TEMP. 
FLAME TEMP. 
RAO. TEMP. 

GAS CONST. 
WALL ABSORPT, 
METAL PCI. 
PART. PCT. 
2ETA 
ZETA-1 
LAMBDA 
DEL-T 
UEL-X/L 
PRINT TIME 
GUIT TIME 
CASE 


1 . 0500 0E . 0 0 
1 .20000E *0 1 
4.27014E-02 
6.30000E*02 
9.10000E*00 
S.00000E.01 
4.33000E-02 

S.40000E *02 
6 , 1 1863E *03 
3.S9164E.03 
3.00000E-03 
9.00000E-01 
2.00000E-01 
0. 

H. 88000E-01 

0 . 

I .B8600E .00 
2.50000E-01 
2.00000E-02 

I , 00000£.00 
S.OOOOOE.OO 
1 



MAR A 1970 


TICE 

0 


OISTANCE 

TOTAL 
REG RATE 

CONVECTIVE 
REG RATE 

OEL-R 

TUBE AREA 

FLOA AREA 

FMISSIVITV 

(X/L) 

(IN/SEC) 

( IN/SEC 1 

(IN) 

(IN**2) 

( I N**2 ) 


1.00000E-01 

7.91429E-03 

1 . 36723E-02 

0. 

B.65901E-01 

7.57831E-01 

4.983948-01 

2.00000E-01 

7 < 31592E-03 

1 .313538-02 

0. 

8.65901E-01 

6.82387E-01 

4.983948-01 

3.00000E-01 

6*9981 IE-03 

1.32913E-02 

0. 

8,6590 1E-C 1 

6. 17856E-0 1 

4 ,983948-0 1 

4 • OOOOOE-O 1 

6.78799E-03 

1.37389E-02 

0. 

8.6590IE-01 

5 .604678-0 1 

4.983948-01 

5 • OOOOUE-O 1 

6.61644E-03 

1 . 369308-02 

0. 

B.65901E-01 

5.40409E-01 

4.98394E-01 

6.00000E-01 

6i47l72E-03 

1.33573E-02 

0. 

B.65901E-01 

5.404098-01 

4,983448-01 

7.00000E-01 

6.353S7E-03 

1.30982E-02 

0. 

8.65901E-01 

5.40409E-01 

4.983048-01 

8.00000E-01 

6.25408E-03 

1.28926E-02 

0. 

B.65901E-01 

5.40409E-01 

4,983948-01 

9.000008-01 

6. 16832E-03 

1 .272648-02 

0. 

B.65901E-OI 

5.40409E-01 

4.983948-01 

1.000008*00 

6;09307E-03 

1.2S901E-02 

0. 

8,6590 1E-0 1 

5. 404098-01 

4.983948-01 

AVG. OEL-R 
(IN) 

0. 

AT. LOSS 
<GM> 

0. 

GAS FLOA 

(LB/SEC) 

5.072088-02 

TOTAL FLO* 
(LB/SEC) 
5.52229E-62 

G-0 

(LB/IN2-SEC) 

5.00057E-02 




tike 

1 


distance 

(X/L) 

1 .00000E-01 
2.00000E-01 
3.00000E-01 
4.00000E-01 
S.00000E-01 
6.00000E-01 
7.00000E-01 
0.000008-01 
9.00000E-01 
1.00000E*00 

AV6. OtL-R 
(IN) 

7.04133E-03 


TOTAL 
REG RATE 
(IN/SEC) 

7.89751E-03 
7 < 30639E- 03 
6.98960E-03 
6.77832E-03 
6.60929E-03 
6.46648E-03 
6.34976E-03 
6.2513TE-03 
6.16648E-03 
6.09193E-03 

AT. LOSS 
(GM> 

5.4349lE*00 


CONVECTIVE 
REG RATE 
(IN/SEC) 

1.33308E-02 

1.28209E-02 

1.29764E-02 

1.3A122E-02 

I.3A291E-02 

I.31072E-02 

1.28590E-02 

1.26624E-02 

1.2S037E-02 

1.23738E-02 

GAS FLOW 
(LB/SEC) 
5.07943E-02 


OEL-R 

(IN) 

7.90782E-03 
7 , 3 1220E-03 
6.99478E-03 
6.78423E-03 
6.61367E-03 
6.46967E-03 
6.3520?E*o3 
6.25299 E-o 3 
6.167578-03 
6.09258E-f>3 

TOTAL FLO* 
(L8/SEC) 
5.S34I0E-S2 


TUBE AREA 
(IN»»2) 

8.92182E-01 

8.90189E-01 

8.8912BE-01 

8.8B424E-01 

8.87854E-01 

8.87374E-01 

8.86981E-01 

B.86650E-01 

8.86365E-01 

8.86115E-01 

G-0 

(LB/IN2-SEC) 

4.86909E-02 


flow area 

I IN*»?) 

7.82113E-01 

7.03482E-01 

6.368888-01 

5.77905E-01 

5.S4U0E-01 

5.S3810E-01 

5.53565E-01 

S.533S8E-01 

5.53180E-01 

5.53024E-01 


FMISSIVITV 


4.98394F-01 

4.983948-01 

4.983948-01 

4.983948-01 

4.983948-01 

4.9B394F-01 

4.983948-01 

4.983948-01 

4.98394E-nl 

4.983448-01 


TIKE 

2 


DISTANCE 

(X/L) 

1.00000E-01 
2.00000E-01 
3.00000E-01 
4 *000008*0 1 
6.00000E-01 
6.00000E-01 
7.00000E-01 
8.00000E-01 
9.00000E-01 
1.000008*00 

AVG. OEL-R 
(IN) 

1 .405948-02 


TOTAL 
REG RATE 
(IN/SEC) 

7.882O0E-O3 

7.29792E-03 

6i98210E-03 

6.76970E-03 

6.60260E-03 

6.46161E-03 

6.34626E-03 

6.24891E-03 

6.164848-03 

6;09093E-03 

AT. LOSS 
(GK) 

1 . 09239E *01 


CONVECTIVE 
REG RATE 
(IN/SEC) 

1 .300378-02 
1.2S190E-02 
1.26739E-02 
1 . 30985E-02 
1.31740E-02 
1 . 286528-02 
1 .262748-02 
1.24393E-02 
1.22878E-02 
1.21641E-02 

GAS FLOW 

(LB/SEC) 

5 . 08687E-02 


OEL-R 

(IN) 

1.57997E-02 

1.461568-02 

1.39817E-02 

1.35594E-02 

1.32206E-02 

1.29345E-02 

1.27007E-02 

1.250368-02 

1.23336E-02 

1.21843E-02 

TOTAL FLO* 
(LB/SEC) 
5.S460SE-62 


TUBE AREA 
(IN*«2) 

9.18803E-01 
9.14784E-01 
9. 12636E-0 1 
9,1 1 206E-0 1 
9. 10060E-01 
9.09093E-01 
9. 08303E-0 1 
9.07637E-01 
9.0T063E-01 
9.06559E-01 

G-0 

(LB/IN2-SEC) 

4.74313E-02 


FLOW AREA 
(IN**2> 

8.06732E-01 
7.24877E-01 
6.S6194E-01 
5.955948-01 
5.67969E-01 
5.67365E-01 
5.66872E-01 
5.664568-01 
5.66098E-0 1 
5.65784E-01 


FMISSIVITV 


4.98272E-01 
4 .98272F-01 
4.98272E-01 
4.982728-01 
4 . 98272E-0 1 
4.98272F-01 
4.98272E-C1 
4. 982728-01 
4.98272F-01 
4.98272E-01 
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TIKE 

3 


DISTANCE 

(X/L) 

1 .OOOOOE-Ol 
2.00000E-01 
3 • 0 OOOOE-O 1 
A.00000E-01 
S.OOOOOE-OI 
6. OOOOOE-Ol 
7. OOOOOE-Ol 
8. OOOOOE-Ol 
9. OOOOOE-Ol 
1.00000E»00 

AVG. OEL-R 
I IN) 

2.10561E-02 


TOTAL 
REG RATE 
(IN/SEC) 

7.86943E-03 

7.29191E-03 

6.97695E-03 

6.76338E-03 

6.59762E-03 

6.4S837E-03 

B.34430E-03 

6.24794E-03 

6.16464E-03 

6.09134E-03 

NT. LOSS 
(GM) 

1.64678E«01 


CONVECT IVE 

reg rate 

(IN/SEC) 

1.26901E-02 
1.22289E-02 
1.23833E-02 
1.27972E-02 
1.29274E-02 
1.26311E-02 
1.24032E-02 
1 ■ 22232E-02 
1.20786E-02 
1.I9608E-02 

GAS FLOW 
(LB/SEC) 
5.09454E-02 


DEL-R 

(IN) 

2.36769E-0Z 
2. 191 1 1E-o2 
2.09617E-o 2 
2.03266E-fl2 
1.98213E-02 
1.93948E-Q2 
1.90461E-02 
1.87S21E-02 
1.84983 E*o 2 
1.82753E-ft2 

total flon 

(LB/SEC) 

5,55838£-n2 


TUBE AREA 
t IN**?) 

9.45764E-01 

9.39686E-01 

9.36427E-0I 

9.34250E-01 

9.32S19E-01 

9.31060E-01 

9.29867E-01 

9.28862E-01 

9.?7996E-0l 

9.27234E-01 

6-0 

(LB/IN2-SECI 

4.62236E-02 


FLOW AREA 
( IN**?) 

R. 31689E-0I 
7 • 46575F-0 1 
6.75776F-01 
6.13535E-01 
5.81985E-01 
5.81074E-01 

S. 80330E-0I 
S.79703E-01 
5.79162E-01 
5.78687E-01 


FMISSIVITY 


4.98272F-01 
4.98272E-P1 
*.9827?E-01 
4.98272F-C1 
4.98272F-01 
4.98272F-01 
4.9827?F-C1 
4.9827?F-0 1 
4.9827?E-01 
4.98272E-M 


TINE 


DISTANCE 

(X/L) 

1 • 0 OOOOE-0 1 
2.00000E-01 
3.00000E-01 
4 • 0 OOOOE-0 1 
5.00000E-01 
6. OOOOOE-O 1 
7 • OOOOOE-O 1 
8.00000E-01 
9. OOOOOE-O 1 
1 . OOOOOE-OO 

AVG. OEL-R 
(IN) 

2.80333E-02 


TOTAL 
REG RATE 
(IN/SEC) 

7.84975E-03 
7.27890E-03 
6.96521E-03 
6.75I02E-03 
6.S863SE-03 
6.44872E-03 
6.33586E-03 
6.24042E-03 
b'. 15785E-03 
6.08514E-03 

NT. LOSS 
(GUI 

2.20674£*01 


convective 

REG HATE 
(IN/SEC) 

1.23890E-02 
1.1 9497E— 02 
1.2I032E-02 
1.25068E-02 
1.26879E-02 
1.24033E-02 
1.21847E-02 
1 .20125E-02 
1 . 18742E-02 
1.1 7620E-02 

GAS FLON 
(LB/SEC) 
S.l 0 1S3E-02 


OEL-R 

(IN) 

3. 15*21E-o2 
2.92012E-02 
2.79372E-02 
2.70881E-02 
2.64173E-I>2 
2.58S22E-02 
2.53<>00E-02 
2.49999E-02 
2.46631E-02 
2.4367JE-02 

TOTAL FLON 
(LB/SEC) 
5.56961E-o2 


TUBE AREA 
(IN**2) 

9.73073E-0I 
9.64905E-01 
9.60508E-01 
9.S7561E-01 
9.55235E-01 
9.53278E-0 1 
9.51679E-01 
9.50331E-01 
9.49167E-0 1 
9,*8145E-01 

G-0 

(LB/IN2-SEC) 

4.50646E-02 


FLON AREA 
(IN**?) 

R.56991E-01 
7.685B1F-01 
6.9563BE-01 
6.31733E-01 
5.961 62E-01 
5.9*9*lE-0l 
5.93943E-01 
5.93101F-01 
5.92375E-01 
5.91737E-01 


FMISSIVITY 


4.96936E-01 

4.96936E-01 

4.96936E-01 

4.96936E-01 

4.96936E-P1 

4.96936E-01 

4.96936F-01 

4.96936E-01 

4.96936E-01 

4.96936E-01 




TIME 

5 


OISTANCE 

(X/L) 

l • OOOOOE-O 1 
2, OOOOOE-O 1 
3 .OOOOOE-O 1 
4.00000E-0! 
5.00000E-01 
6.00000E-01 
7 . OOOOOE-Ol 
8. 00000E-01 
9. OOOOOE-Ol 

i.oooooe*oo 

AVG. OEL-R 
(IN) 

3.49873E-02 


TOTAL 
REG RATE 
(IN/SEC) 

7.84082E-03 
7.27590E-03 
6.96294E-03 
6.74759E-03 
6 . 58324E-03 
6.44719E-03 
6 . 33549E-03 
6.24094E-03 
6. 15905E-03 
6.08687E-03 

NT. LOSS 
(GM) 

2.77192E*01 


CONVECTIVE 
REG RATE 
(IN/SEC) 

1 .21002E-02 
1 . 16815E-02 
1.1 8343E-02 
1 .22284E-02 
1.24572E-02 
1 .2 1839E-02 
1 , 19743E-02 
1.18095E-02 
1 . 16775E-02 
1.1S706E-02 

GAS FLON 

(LB/SEC) 

5.10949E-02 


DEL-R 

(IN) 

3.93B83E-02 

3.64789E-02 

3.4901SE-02 

3.38377E-02 

3.30024F-02 

3.23003E-02 

3.17256E-02 

3.12405E-02 

3.08213E-02 

3.04528E-02 

TOTAL FLON 
(LB/SEC) 
5.58239E-S2 


TUBE AREA 
(IN**?) 

l.00070E*00 
9.9041 3E-0 1 
9.84856E-01 
9.8111 7E-0 1 
9.781 87E-0 1 
9.75726E-01 
9.73715E-01 
9.72019E-C1 
9.70555E-01 
9.69268E-01 

G-0 

(LB/IN2-SEC) 

4.39523E-02 


FLON AREA 
(IN**?) 

8.82613F-01 
7.90874F-01 
7 . 15762E-0 1 
6.501 73E-01 
6.1 0486E-0 1 
6 . 08951E-0 1 
6.07696E-01 
6 • 06637E-0 1 
6.05723E-01 
6.04920E-01 


FMISSIVITY 


4.96936F-P1 

4.96936E-P1 

4.96936E-01 

4.96936F-P) 

4.96936E-01 

4.96936E-01 

4.96936E-C1 

4.96936E-01 

4.9693hF-OI 

4.96936E-C1 


36 



REFERENCES 


1. Marxman, G. A. , Wooldridge, C. E. , and Muzzy, R. J. , "Fundamentals 
of Hybrid Boundary Layer Combustion," Heterogeneous Combustion , 
Progress in Astronautics and Aeronautics , Vol. 15, Academic Press, 

New York, 1964, p. 485. 

2. Marxman, G. A., and Gilbert, M. , "Turbulent Boundary Layer Combustion 
in the Hybrid Rocket," Ninth Symposium (International) on Combustion , 
Academic Press, New York, 1963, p. 371. 

3. Barrere, M. , and Moutet, A., "La Propulsion par Fusees Hybrides," 
presented at the International Astronautical Congress, Paris, Sep- 
tember 1963. 

4. Green, L. , Jr. , "introductory Considerations of Hybrid Rocket Com- 
bustion," Heterogeneous Combustion, Progress in Astronautics and 
Aeronautics , Vol. 15, Academic Press, New York, 1964, p. 451. 

5. Marxman, G. A., and Wooldridge, C. E. , "Research on the Combustion 
Mechanism of Hybrid Rockets," presented at the Combustion and Pro- 
pulsion Panel, AGARD, NATO, San Diego, California, April 22-24, 

1965. 

6. Marxman, G. A. , "Combustion in the Turbulent Boundary Layer on a 
Vaporizing Surface," Tenth Symposium (International) on Combustion , 
The Combustion Institute, Pittsburgh, 1965, p. 1337. 

7. Wooldridge, C. E. , Marxman, G. A., and Kier, R. J. , "investigation 
of Combustion Instability in Hybrid Rockets", Final Report, NASA 
Contract NAS 1-7310, Stanford Research Institute, September 1969. 

8. Smoot, L. D. , and Price, C. F. , "Regression Rate Mechanism of Non- 
metalized Hybrid Fuel Systems," AIAA J. 3, 1408-1413 (1965). 

9. , "Regression Rates of Metalized Hybrid Fuel Systems," 

AIAA J. 4, 910-915 (1966). 

10. , "Pressure Dependence of Hybrid Fuel Regression Rates," 

AIAA J. 5, 102-106 (1967). 


37 



